The paradox of soft singularity crossing and its resolution 
by distributional cosmological quantitities 
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A cosmological model of a fiat Friedmann universe filled with a mixture of anti-Chaplygin gas 
and dust-like matter exhibits a future soft singularity, where the pressure of the anti-Chaplygin gas 
diverges (while its energy density is finite). Despite infinite tidal forces the geodesies pass through 
the singularity. Due to the dust component, the Hubble parameter has a non-zero value at the en- 
counter with the singularity, therefore the dust implies further expansion. With continued expansion 
however, the energy density and the pressure of the anti-Chaplygin gas would become ill-defined, 
hence from the point of view of the anti-Chaplygin gas only a contraction is allowed. Paradoxi- 
cally, the universe in this cosmological model would have to expand and contract simultaneously. 
This obviosly could not happen. We solve the paradox by redefining the anti-Chaplygin gas in a 
distributional sense. Then a contraction could follow the expansion phase at the singularity at the 
price of a jump in the Hubble parameter. Although such an abrupt change is not common in any 
cosmological evolution, we explicitly show that the set of Friedmann, Raychaudhuri and continuity 
equations are all obeyed both at the singularity and in its vicinity. We also prove that the Israel 
junction conditions are obeyed through the singular spatial hypersurface. In particular we enounce 
and prove a more general form of the Lanczos equation. 



PACS numbers: 98.80. Cq, 98.80. Jk, 98.80.Es, 95.36.-l-x 



I. INTRODUCTION 



The problem of cosmological singularities has been at- 
tracting the attention of theoreticians working in gravity 
and cosmology since the early fifties In the sixties 

general theorems about the conditions for the appear- 
ance of singularities were proven [J, and the oscillatory 
regime of approaching the singularity Q , the Mixmaster 
universe 0] was discovered. Basically, until the end of 
nineties almost all discussions about singularities were 
devoted to the Big Bang and Big Crunch singularities, 
which are characterized by a vanishing cosmological ra- 
dius. 

However, kinematical investigations of Friedmann cos- 
mologies have raised the possibility of sudden future sin- 
gularity occurrence characterized by a diverging 
a whereas both the scale factor a and a are finite. Then, 
the Hubble parameter H = a/a and the energy density 
p are also finite, while the first derivative of the Hub- 
ble parameter and the pressure p diverge. Until recent 
years, however, sudden future singularities attracted only 
a limited interest among researchers. The interest grew 
due to two reasons. The recent discovery of the cosmic 
acceleration [l^ has stimulated the elaboration of dark 
energy models, responsible for such a phenomenon (see 
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e.g. for review [20] )• Remarkably in some of these mod- 
els the sudden singularities arise quite naturally. Another 
source of the interest to sudden singularities is the devel- 
opment of brane models [lol . [ill . [l8j , where singularities 
of this kind could arise naturally (sometimes these sin- 
gularities, arising in brane- world models, are called "qui- 
escent" [3). 

In the investigations devoted to sudden singularities 
one can distinguish three main topics. The first of them 
deals with the question of the compatibility of the mod- 
els possessing soft singularities with observational data 
[isl I21L [2^ 1 . The second direction is connected with the 
study of quantum effects [III [13 l23l - [25| . Here one can 
see two subdirections: the study of quantum corrections 
to the effective Friedmann equation, which can elimi- 
nate classical singularities or at least, change their form 
[lol [TtI . 1 2 3] ; and the study of solutions of the Wheeler- 
DeWitt equation for the quantum state of the universe 
in the presence of sudden singularities [23 . [25j . The third 
direction is connected with the possibility of the sudden 
singularity crossing in classical cosmology (25l - [29j . The 
present paper is devoted exactly to this topic. 

A particular feature of the sudden future singularities 
is their softness [1^. As the ChristofFel symbols depend 
only on the first derivative of the scale factor, they are 
regular at these singularities. Hence, the geodesies are 
well behaved and they can cross the singularity [IB] • One 
can argue that the particles crossing the singularity will 
generate the geometry of the spacetime, providing in such 
a way a "soft rebirth" of the universe after the singularity 
crossing [2^ . Note that the possibility of crossing of some 
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kind of cosmological singularities was noticed already in 
the early paper by Tipler [^fl] • A close idea of integrable 
singularities in black holes, which can give origin to a 
cosmogenesis was recently put forward in [3l'|. 

As a starting point we consider an interesting exam- 
ple of a sudden future singularity - the Big Brake which 
was discovered in Ref. [33| while studying a particular 
tachyon cosmological model. The particularity of the Big 
Brake singularity consists in the fact that the time deriva- 
tive of the scale factor is not only finite, but exactly equal 
to zero. That makes the analysis of the behavior in the 
vicinity of singularity especially convenient. In particu- 
lar, in Ref. [23l it was shown that the predictions of the 
future of the universe in this model [32l | are compatible 
with the supernovae type la data, while in Refs. jlJ, HI] 
some quantum cosmological questions were studied in the 
presence of the Big Brake singularity. 

The simplest cosmological model allowing a Big Brake 
singularity was also introduced in Ref. [32J. This model 
is based on the perfect fluid, dubbed "anti-Chaplygin 
gas" . This fluid is characterized by the equation of state 



where A is a positive constant. Such and equation of 
state arises, for example, in the theory of wiggly strings 
[33} . In paper [s?] a fluid obeying the equation of state 
^ was called "anti-Chaplygin gas" in analogy with the 
Chaplygin gas [34] which has the equation of state p = 
—A/p and has acquired some popularity as a candidate 
for a unified theory of dark energy and dark matter [s^ . 

An explicit example of the crossing of the Big Brake 
singularity was described in detail in paper [2^ , were the 
tachyon model [13, [s^l was investigated. In this model 
the tachyon field passes through the singularity, continu- 
ing its evolution with a recoUapse towards a Big Crunch. 
In a simpler model, based on the anti-Chaplygin gas, such 
a crossing is even easier to understand. 

The next natural step in the analysis of the soft sin- 
gularities seems to be obvious. One can consider a soft 
singularity of more general type than the big brake by 
adding to the tachyon matter or to the anti-Chaplygin 
gas some dustlike matter. However, in this case the 
traversability of the singularity seems to be obstructed. 
The main reason for this is that while the energy den- 
sity of the tachyonic field (or of the anti-Chaplygin gas) 
vanishes at the singularity, the energy density of the mat- 
ter component does not, leaving the Hubble parameter 
at the singularity with a finite value. Then some kind 
of the paradox arises: if the universe continues its ex- 
pansion, and if the equation of state of the component 
of matter, responsible for the appearance of the soft sin- 
gularity (in the simplest case, the anti-Chaplygin gas) is 
unchanged, then the expression for the energy density of 
this component becomes imaginary, which is unaccept- 
able. The situation looks rather strange: indeed, the 
model, including dust should be in some sense more reg- 
ular, than a single exotic fluid, the anti-Chaplygin gas. 



Thus, if the model based on the pure anti-Chaplygin gas 
has a traversable Big Brake singularity, than the more 
general singularity arising in the model based on the mix- 
ture of the anti-Chaplygin gas and dust should also be 
traversable. 

Related to that, it was recently shown that general soft 
singularities arising in the Friedmann model, filled with 
the scalar field with a negative potential, inversely pro- 
portional to this field are traversable. So, what could be 
wrong with the simple two-fiuid model? One can see that 
what we face is some sort of a clash between the equa- 
tion of state of one of these fluids and the dynamics (the 
Friedmann and Raychaudhuri equations) and energy con- 
servation equations. In this paper we shall try to resolve 
this paradox, insisting on the preservation of the equa- 
tion of state of the anti-Chaplygin gas. The price which 
one has to pay for it is the obligatory use of the general- 
ized functions for some cosmological quantities. Namely, 
the anti-Chaplygin gas remains physical if rather a recol- 
lapse follows, but then the Hubble parameter would have 
a sharp jump, obstructing the validity of the Raychaud- 
huri equation (the second Friedmann equation) in the 
usual sense of functions. Thus, apparently, the evolution 
cannot be continued through the soft singularity, unless 
treating the cosmological quantities as distributions. We 
claim that such a generalization is mathematically rig- 
orous, moreover, the introduction of distributions is not 
so drastic, as it looks at the first sight, as the pressure 
of the anti-Chaplygin gas diverges anyhow at the soft 
singularity (as it so does for the tachyon held). Then 
in the Conclusion we shall dwell on the possible physi- 
cal sense of the proposed constructions and its possible 
alternatives. 

The plan of the paper is as follows. In Section [H] 
we discuss generic Friedmann space-times, which admit 
H — —00 type singularities, while the Hubble parameter 
H remains finite. Such singularities are related to corre- 
sponding divergencies in the pressure of the perfect fiuid 
filling the Friedmann universe (while its energy density 
stays finite). We investigate the kinematics, the geodesic 
equations, the geodesic deviation equations in the vicin- 
ity of these singularities and also prove that these singu- 
larities are weak. 

In Section IIIII we discuss a mixture of the anti- 
Chaplygin gas and dust in a flat Friedmann universe and 
explain the essence of the paradox. We explicitly derive 
the behavior of the energy density and pressure in the 
vicinity of the soft singularity and we solve the geodesic 
equations in this region. While the singularity turns to 
be traversable by the geodesies, the explicit solution also 
shows that the Raychaudhuri equation is violated at the 
singularity. 

In Section IIVI we add generalized distributional con- 
tributions to both the pressure and energy density, such 
that (a) the equation of state of the anti-Chaplygin gas 
still holds and (b) the singularity becomes traversable. 
We also perform checks of the Friedmann, Raychaudhuri 
and continuity equations, which all hold valid across the 
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singularity in a distributional sense. In the process we 
employ a number of Propositions on distributions pre- 
sented and proved in Appendix A. For the convenience 
of the reader we present a related semi-heuristic discus- 
sion of two known distributional identities in Appendix 
B. We stress that the distributional modifications of the 
energy density and pressure do not modify the cosmolog- 
ical evolution, but they make possible the soft singularity 
crossing. 

In Section |V] we revisit the junction conditions along 
a spacelike hypersurface in a flat Fricdmann universe. 
The future soft singularity represents such a spatial hy- 
persurface, along which the energy-momentum tensor di- 
verges. Extending the space-time through this hypersur- 
face is possible by obeying both Israel junction condi- 
tions. While the first condition, requiring the continu- 
ity of the induced metric is easy to satisfy (the metric 
stays regular at the soft singularity), the second con- 
dition relates the jump in the extrinsic curvature to 
the distributional part of the energy-momentum tensor 
through the Lanczos equation. We will show that in flat 
Friedmann space-times the Lanczos equation holds for 
a more general class of distributional energy-momentum 
tensors. With this we give a second proof that the gen- 
eralized distributional energy-momentum tensor assures 
the traversability of the soft singularity. In the process 
we employ a simple form of the Lanczos equation valid 
in flat Friedmann universes, derived in Appendix [Cl 

We summarize our results and give some further out- 
look in the Concluding Remarks. 

We chose c = 1 and SttG/S = 1. A subscript S denotes 
the value of the respective quantity at the soft singularity. 



II. PRESSURE SINGULARITIES IN FLAT 
FRIEDMANN UNIVERSES 

The line element squared of a flat Friedmann universe 
can be written as 



ds^ = -dt^ + a^ it)Y,{dx°') 



(2) 



where {a = 1,2,3) are Cartesian coordinates. The 
evolution of the Friedmann universe is governed by the 
Raychaudhuri (second Friedmann) equation 

Up + p) ^ (3) 



H 



and by the continuity equation for the fluid 

p + 3H{p + p)=0. (4) 

Here the dot denotes the derivative with respect to cos- 
mological time t. A first integral of this system is given 
by the first Friedmann equation 



H^=p 



(5) 



It is easy to see that the Raychaudhuri equation can 
be obtained from the first Friedmann and the continu- 
ity equations. 



A. Kinematics in the vicinity of sudden 
singularities 

Sudden singularities are characterized by finite H and 
— > — oo (finite d and d — > —oo) at some finite scale 
factor a. The energy density of the fluid is finite but its 
pressure diverges at this type of singularity, therefore the 
term "pressure singularity" is also in use. Then, we would 
like to emphasize the fact that there is an essential dif- 
ference between the sudden singularities with H — and 
with H > 0. As has been already mentioned in the Intro- 
duction, in the first case, which is called Big Brake, the 
universe begin contracting and running towards the Big 
Crunch singularity. Exactly this occurs in models based 
on tachyon field with a particular potential [2^ . [29l l32j 
or in the anti-Chaplygin gas models. In the case of the 
model based on the mixture of one of this fluids and dust, 
we encounter the second situation when the value of the 
Hubble constant is positive at the moment of encounter 
with the sudden singularity. That means that after cross- 
ing the singularity the universe should continue its expan- 
sion, but the anti-Chaplygin gas becomes ill-defined, as 
it will be shown in detail in Section Hill devoted to the 
model based on mixture of the anti-Chaplygin gas and 
dust. 

One possible way of overcoming this obstacle is to al- 
low the jump in the sign of the Hubble parameter, which 
as was mentioned in the Introduction leaves valid the 
first Friedmann equation, the continuity equation and the 
equation of state, while making invalid the Raychaudhuri 
equation. This last obstacle can be cured by the accep- 
tance the distributional Dirac (5-function type contribu- 
tions into the pressure and the energy density, which is 
described in detail in the section IV. 



B. Geodesies in the vicinity of sudden singularities 

The geodesic equations in flat Friedmann space-time 
are 



2-- 



dX a dX dX 



, 



d'^t / dx" ' , 

ad 2^ I — I =0 



dX' 



dX 



(6) 



(7) 



where A is an affine parameter. Integrating these equa- 
tions yields 



dx" _ P" 
~dX~~i?' 



dtV 

= €+ — 



dX 



(8) 



(9) 



with P", e integration constants and P^ = J^a {P")^- 
The quantity e is fixed by the length of the tangent vector 
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u°' of the geodesic as e = —UaU"", i.e. one for timelike and 
zero for lightlike orbits. In a comoving system P" — 
and t = A is affine parameter. 

Eqs. ([8]) and (O are singular only for vanishing scale 
factor (see also Ref. [1^). Therefore, the existence of a 
solution t (A), x" (A) of Eqs. ^ and ^ is assured by the 
Cauchy-Peano theorem for any nonzero a (including the 
soft singularity). Thus the functions t (A) and x°' (A), i.e. 
the geodesies can be continued through the singularity 
occurring at finite scale factor. Only derivatives of higher 
order than two of t (A) and (A) are singular (as they 
contain a), however these do not appear in the geodesic 
equations. Pointlike particles moving on geodesies do 
not experience any singularity. Thus, as we argued in 
the preceding paper [29] one is not obliged to consider 
such a singularity as a final state of the universe. Indeed, 
passing through this singularity the matter recreates also 
the spacetime in a unique way, at least for such simple 
models, as those based on Friedmann metrics. 



C. Deviation equation in the vicinity of sudden 
singularities 

The 3-spaces with t =const have vanishing Riemann 
curvature. However, the 4-dimensional Riemann curva- 
ture tensor has the nonvanishing components: 

212 = 313 = R^ 323 = (1^) 

and the corresponding components arising from symme- 
try. Here a,/3 = 1,2,3. Remarkably, all components 
which diverge at the singularity are of the type Rtata 
[29j . Therefore, the singularity arises in the mixed spatio- 
temporal components. 

The geodesic deviation equation along the integral 
curves of u = d / dt (which are geodesies with affine pa- 
rameter t) is 

it'^ = -i?^,,r;Vw^ (11) 

where rf' is the deviation vector separating neighboring 
geodesies, chosen to satisfy ifub — 0. For a Friedmann 
universe it becomes 

u- = -R\,^rf ^h, (12) 

which at the singularity diverges as — oo. Therefore, 
when approaching the singularity, the tidal forces mani- 
fest themselves as an infinite braking force stopping the 
further increase of the separation of geodesies, but not 
the evolution along the geodesies. With a < in the 
vicinity of the singularity, once the geodesies have passed 
through, they will approach each other. Therefore a con- 
traction phase will follow: everything that has reached 
the singularity will bounce back. 



D. The type of the singularity 

In this subsection we shall present the classification 
of singularities, based on the point of view of finite size 
objects, which approach these singularities. In princi- 
ple, finite size objects could be destroyed while passing 
through the singularity due to the occurring infinite tidal 
forces. A strong curvature singularity is defined by the 
requirement that an extended finite object is crushed to 
zero volume by tidal forces. We give below Tipler's [s^l 
and Krolak's ^] definitions of strong curvature singular- 
ities together with the relative necessary and sufficient 
conditions. An alternative definition of the softness of 
a singularity, based on a Raychaudhuri averaging, was 
developed by Dabrowski js^] . 

According to Tipler's definition if every volume ele- 
ment defined by three linearly independent, vorticity- 
free, geodesic deviation vectors along every causal 
geodesic through a point p vanishes, a strong curvature 
sing ularity is encountered at the respective point p [3d| , 
|26| . The necessary and sufficient condition for a causal 
geodesic to run into a strong singularity at As (A is affine 
parameter of the curve) ^38,] is that the double integral 

/ dX' / dX" \R\ji,u''u''\ (13) 
Jo Jo 

diverges as A As . A similar condition is valid for light- 
like geodesies, with R^^^f^u'^u'' replacing Rabu'^u^ in the 
double integral. 

Krolak's definition is less restrictive. A future-endless, 
future-incomplete null (timelike) geodesic 7 is said to ter- 
minate in the future at a strong curvature singularity 
if, for each point p G 7, the expansion of every future- 
directed congruence of null (timelike) geodesies emanat- 
ing from p and containing 7 becomes negative somewhere 
on 7 [3^ , [3^ . The necessary and sufficient condition for 
a causal geodesic to run into a strong singularity at As 
[38t is that the integral 

[\x' \R\^,u^u''\ (14) 
Jo 

diverges as A As. Again, a similar condition is valid 
for lightlike geodesies, with replacing Rabu'^u^ 

in the integral. 

In flat Friedmann space-time the comoving observers 
move on geodesies having four velocity u = d/dt, where 
t is affine parameter. The nonvanishing components of 
Riemann tensor are given by Eq. (|10p . Since H is finite 
along the geodesies, neither of the integrals (|13p and 
diverge at the singularity. The singularity is weak (soft) 
according to both Tipler's and Krolak's definitions. That 
means although the tidal forces become infinite, the finite 
objects are not necessarily crushed when reaching the 
singularity (see also f26j). 
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III. THE PARADOX OF THE SOFT 
SINGULARITY CROSSING IN THE 
COSMOLOGICAL MODEL BASED ON THE 
ANTI-CHAPLYGIN GAS AND DUST UNIVERSE 



is the value of the Hubble parameter at t = ts- Corre- 
spondingly the leading terms of the energy densities of 
the anti-Chaplygin gas and dust, also of the pressure of 
the anti-Chaplygin gas are 



We discuss an universe filled with two components. 
One is the anti-Chaplygin gas with the equation of state 
((H) and other is the pressureless dust. 

The solution of the continuity equation for the anti- 
Chaplygin gas gives the following dependence of its en- 
ergy density on the scale factor: 



PACh 



4-A 



(15) 



where i? is a positive constant, determining the initial 
condition. The energy density of the dust-like matter is 
as usual 



Pm,0 



(16) 



where g is a constant. 

It is clear that when during the expansion of the uni- 
verse, its scale factor approaches the value 



PACh 



PACh 



^6AHs{ts-t) 



I A 
6Hs{ts-t) 



(20) 



(21) 



(22) 



One can see that the expressions p^ . (PT|) and ([^ 
cannot be continued for t > ts due to the emerging neg- 
ative quantities under the square roots. The assumption 
of a sharp transition from expansion to contraction im- 
plies the following changes in Eqs. (IT51) - (1^^ : 



a(t) = as 



Pin, a 

as 



ts-t\- 



l2Aal, 



|3/2 



(23) 



B 



(17) 



P,n^Hl + 3Hl\ts-t\ , 



the energy density of the anti-Chaplygin gas vanishes, 
and its pressure grows to infinity. That means that the 
deceleration also becomes infinite. However, the energy 
density of dust remains finite, hence the same is true also 
for the Hubble parameter. It is here that the paradox 
arises: if the universe continues to expand, the expression 
under the sign of the square root in Eq. ([T5l) becomes 
negative and the energy density of the anti-Chaplygin 
gas becomes ill-defined. A way out of this situation is 
only by assuming that at this moment the Hubble pa- 
rameter changes its sign, while keeping its absolute value 
(such that the energy density will not have a jump, as im- 
plied by the Friedmann equation). This possibility will 
be studied in detail in the following subsections. 



A. Evolutions in the vicinity of the singularity 



p^ch = V(i^Hs\ts-t\ 



A 



PACh 



GHs\ts~t\ 



(24) 



(25) 



(26) 



The quantities are well-defined and continuous 

at the moment of the singularity crossing. The expres- 
sion for the pressure (|26|) is divergent, but this diver- 
gence is integrable and this is sufficient for our purposes. 
These new expressions satisfy the Friedmann equation, 
the continuity equations and the equation of state for 
the anti-Chaplygin gas. However, the time derivatives of 
these quantities are not continuous and it is the reason 
of the failure of the Raychaudhuri equation. We shall 
analyze this problem in the following section, but before 
we discuss the geodesies in the vicinity of the singularity. 



Let us substitute the expressions ([T5|) and ([T6| into the 
first Friedmann equation. We shall find its solution for 
the universe approaching to the soft singularity point at 
the moment ts (the latter cannot be found analytically, 
but its value is not important for our analysis): 



a{t) = as 
where 



Pm,Q 

as 



{ts -t)- 



1 2Aal , 



3Hs 



Hs 



Pm.O 



tf'^ , (18) 



(19) 



B. Singularity crossing geodesies 

We can integrate explicitly the geodesies equations ^ 
and ([9]) in the vicinity of singularity, using the expression 
(1^51) for the cosmological factor, also taken in the vicinity 
of singularity. Choosing the affine parameter in such a 
way that the point A = corresponds to the singularity 
crossing we obtain up to the second order in A terms 



t^ts 




2a| 



(27) 
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, P' P'^Hs ,,,,2 
e H 2 2 — sgn[A)A 



(28) 



One can see from Eqs. (P7)) and (|28l) that not only the 
time and spatial coordinates of the geodesies are contin- 
uous at the soft singularity crossing, but also their first 
derivatives with respect to the afhne parameter A. 



IV. SINGULARITY CROSSING, THE 
RAYCHAUDHURI EQUATION AND 
DISTRIBUTIONS 

Let us discuss the expressions for the Hubble parame- 
ter and its time derivative in the vicinity of the singular- 
ity. Starting from the expression (I23p we obtain 

H{t) = Hssgn{ts-t) 

/ 



^ sgn{ts 



t)^\ts-t\ , (29) 



H = -2Hs5{ts - t) 



I 3A sgnits - 1) 



(30) 



Naturally, the 5-term in H arises because of the jump 
in H , as the expansion of the universe is followed by a 
contraction. To restore the validity of the Raychaudhuri 
equation we shall add a singular (5-terni to the pressure 
of the anti-Chaplygin gas, which will acquire the form 



PACh 



The equation of state ^ of the anti-Chaplygin gas is 
preserved, if we also modify the expression for its energy 
density: 



A 



PACh 



6Hs\ts-t\ 



\HsSits - 1) 



(32) 



The last expression should be understood in the sense 
of the composition of distributions (see Appendix A and 
the references therein). 

In order to prove that pACh and p^^^ represent a self- 
consistent solution of the system of cosmological equa- 
tions, we shall use the following distributional identities: 



[sgn{T)g{\T\)]6{T) = , 



[f{r)+CS{T)] 



r\r) , 



±[fir)+C5{r)]-' = ^rHr) 

ClT ClT 



(33) 
(34) 

(35) 



Here 5 (|t|) is bounded on every finite interval, / (t) > 
and C > is a constant. These identities follow from the 
Propositions 1, 2 and the Corollary enounced and proved 



in the Appendix A. The parameter r stays instead of the 
difference ts — t. 

Because of Eq. the energy density ([5^ behaves as 
a continuous function which vanishes at the singularity. 
The first term in the expression for the pressure (j3ip 
diverges at the singularity. Therefore the addition of a 
Dirac delta term, which is not changing the value of pACh 
at any r 7^ (i.e. t ^ ^5) does not look too drastic and 
might be considered as a some kind of renormalization. 

To prove that Friedmann, Raychaudhuri and continu- 
ity equations are satisfied we must only investigate those 
terms, appearing in the field equations, which contain 
Dirac (5-functions, since without them, these equations 
can be reduced to those we have found in the previ- 
ous section. First, we check the continuity equation for 
the anti-Chaplygin gas. Due to the identities 
the S (r)-terms occurring in p^^/j and Pj^(jf^ could be 
dropped. We keep them however in order to have the 
equation of state explicitly satisfied. Then the 5 (T)-term 
appearing in 2>HpACh vanishes, because the Hubble pa- 
rameter changes sign at the singularity [see Eq. ([55]) ] . 

The 5 (r)-term appearing in p^ch does not affect the 
Friedmann equation due to the identity (j34l) . Finally, the 
(5-term arising in the time derivative of the Hubble pa- 
rameter in the left-hand side of the Raychaudhuri equa- 
tion is compensated by the conveniently chosen (5-term 
in the right-hand side of Eq. (PT|) . 



V. THE JUNCTION CONDITIONS ACROSS 
THE SINGULARITY 

In this section we discuss the singularity crossing in a 
slightly different way, by analyzing the junction condi- 
tions. We have to match two space-time regions across 
the space-like hypersurface t — Q. The junction of two 
space-time regions has to obey the Israel matching condi- 
tions 40] , namely, the induced metric should be continu- 
ous and the extrinsic curvature of the junction hypersur- 
face could possibly have a jump, which is related to the 
distributional energy-momentum tensor on the hypersur- 
face by the Lanczos equation. The scale factor being con- 
tinuous across the singularity the first Israel condition is 
obeyed. We will next prove that the second Israel junc- 
tion condition (the Lanczos equation [4lj, [i^]) is also 
satisfied. 

For this we have to check whether Eq. 



(36) 



(see Appendix [Ct . still implies the Lanczos equation 

Ai?k = 'Ip , (37) 

derived in Appendix [Cl when p ^ p5 (r) — pACh, p — 
Pm + PACh and PACh generalized to a distribution 

P(r) 



PACh^ 



[R{T)^Q{T)5{r)r 



(38) 
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Here uj > 0, R{t) > 0, Q (r) > and P (r) is bounded. 

When Eq. ([5^1) is applied to a test function ip (r), the 
terms containing and p„ give regular contributions 
and the limits of the respective integrals vanish, similarly 
as discussed in Appendix ICl Also, due to Proposition 2 
given in the Appendix A, the integral of the distributional 
term containing p^c'/j becomes 



i?" (r) 



dT 



(39) 



which also vanishes for e 
the contributions 



0. Wc still have to consider 



' (r)] If (r) a^dr 



(40) 



Although the contribution pip (r) to the integrand is 
singular at r = 0, its integral can be conveniently evalu- 
ated by the Residue Theorem. For this we remark, that 
the integrand is an analytically extendible function into 
the complex plane in the vicinity of r = and its residue 
is zero, therefore the integral vanishes. Finally, the con- 
tribution containing p leads by integration and the limit- 
ing process to the right hand side of the Lanczos equation 
(El). 

Therefore we have proven that the space-time regions 
separated by the singular spatial hypersurface, represent- 
ing the pressure singularity, can be joined. In other 
words, the singularity becomes traversable. 



VI. CONCLUDING REMARKS 

It is known that certain models of cosmological fluids 
in Friedmann universes, like the anti-Chaplygin gas or 
the tachyon field with a special potential [33 , evolve into 
a sudden future singularity, which in spite of a diverg- 
ing pressure, is weak. It was argued that singularities of 
this kind could be traversable despite infinite tidal forces 
emerging at the singularity for an infinitesimally short 
time [26]. In Ref. 29] the process of crossing of the 
Big Brake singularity was described in some detail for 
the tachyon model [S^]- (The particularity of the Big 
Brake singularity, consists in the fact that at the cross- 
ing of such a singularity the Hubble variable is not only 
finite, but vanishes.) We also note recent discussions [42| 
on crossing the "traditional" Big Bang and Big Crunch 
singularities. 

In the present paper we considered a simple cosmologi- 
cal model containing a mixture of anti-Chaplygin gas and 
dust. We have shown that the geodesies equations and 
their solutions are still well-defined in this case, however 
the inclusion of dust generates a nonzero value of the 
Hubble parameter at the singularity encounter, gener- 
ating the following paradox. The dust would require a 
continued expansion, which would make the energy den- 
sity and pressure of the anti-Chaplygin gas ill-defined. A 



contraction in turn, would be compatible with the anti- 
Chaplygin gas, nevertheless implying an abrupt change 
of the Hubble parameter from expansion to contraction. 
The jump in the Hubble parameter implies the appear- 
ance of the 6 function in the Raychaudhuri equation 
(which contains H). 

We have cured this situation by redefining the pressure 
and energy density of the anti-Chaplygin gas as distribu- 
tions. As an equivalent interpretation, the pressure can 
be generalized by the addition of a distributional con- 
tribution, while the energy density left unchanged, at 
the price of redefining the equation of state of the anti- 
Chaplygin gas in a distributional sense. Then all cosmo- 
logical equations are satisfied in the same distributional 
sense. We have also shown, that the Israel junction con- 
ditions are obeyed through the singular spatial hypersur- 
face, in particular we have enounced and proved a more 
general form of the Lanczos equation. The results rely 
on two Propositions and a Corollary proven in Appendix 
A. 

The resolution of the paradox at the soft singularity 
crossing by the introduction of distributional quantities 
and equations may look unusual, however distributional 
quantities, localized on hypersurfaces are quite com- 
monly used in general relativity and other gravitational 
theories. Spacetime regions are frequently matched by 
the inclusion of distributional layers; also shock-waves 
can be modeled by Dirac (5-functions. Braneworld mod- 
els [ill, [111 arise due to the orbifold boundary condi- 
tions, the non-smoothness of the 5-dimensional metric 
at the brane (the jump in its extrinsic curvature) be- 
ing directly related to the distributional 3-1-1 standard 
model fields embedded in the 5-dimensional spacetime. 
Besides, metrics allowing distributional curvature were 
considered earlier for studying strings and other distri- 
butional sources in general relativity [i^. The appli- 
cations of the distributional quantities to the study of 
Schwarzschild geometry and point massive par ticles in 
general relativity were used in Refs. [1^ and 47 1 respec- 
tively. 

More generically the connection between singularities 
and the distributional treatment of the physical quanti- 
ties is well-known in quantum field theory. Indeed, the 
appearance of the ultraviolet divergences can be under- 
stood as the result of the indefiniteness of the product 
of distributions and the renormalization procedure could 
be interpreted as a definition of such a product [48j . 

We hope that the investigations presented here may 
turn useful in deriving similar results in connection with 
the traversability of other types of sudden singularities. 

While mathematically self-consistent, the scenario pre- 
sented in this paper may look somewhat counter-intuitive 
from the physical point of view. This is because its es- 
sential ingredient is the abrupt change of the expansion 
into a contraction. However, such a behavior is not more 
counter-intuitive that the absolutely elastic bounce of the 
ball from a rigid wall, as known in classical mechanics. 
Indeed, in the latter case the velocity and the momen- 
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turn of the ball change their direction abruptly. That 
means that an infinite force acts from the wall onto the 
ball during an infinitely small interval of time. The re- 
sult of this action is however integrable and results in a 
finite change of the momentum of the ball. In fact, the 
absolutely elastic bounce is an idealization of a process 
of finite time-span during which inelastic deformations of 
both the ball and the wall are likely. It is reasonable to 
think that something similar occurs also in the two-fiuid 
universe model presented in this paper, which undergoes 
a transition from an expanding to a contracting phase. 
The smoothing of this process should involve some (tem- 
porary) geometrically induced change of the equation of 
state of matter. Note, that such changes are not un- 
common in cosmology. In the tachyon model '32] which 
was starting point of our studies of the Big Brake sin- 
gularities, there was the tachyon -pseudotachyon trans- 
formation driven by the continuity of the cosmological 
evolution. In a cosmological model with the phantom 
field with the cusped potential [i^, the transformations 
between phantom and standard scalar field were consid- 
ered. Thus, one can imagine that the real process of the 
transition from the expansion to contraction induced by 
passing through a soft singularity can imply some tem- 
porary change of the equation of state which makes the 
above processes smoother. We hope to explore such a 
scenario in the future. 



functions having continuous derivatives of all orders and 
compact support. The action of an / 6 2?' distribu- 
tion on test functions is given by (/, Lp) , which in the 
case when / is an ordinary locally summable function 
is nothing but / (r) ip (r) dr. We note that (5„ (r) 
has the compact support: 1/n). We will also use 

the n-th derivative of / G 2?' acts as {df (t) /dT",(p) — 
(-l)"(/(r),d^/dr"). 

For an arbitrary distribution /, the function /„ (r) — 
f * Sn = if {t — x) , Sn (x)) gives a sequence of infinitely 
differentiable functions converging to /. 

Definition 1. ; The commutative product of f and g 
exists and is equal to h on the open interval (a, b) (~oo < 
a < b < oo) if 

lim {fn9n,'P) = {h,^} 
n— >oo 

for any Lp G T> with support contained in the interval 
ia,b) fM]'- 

Proposition 1. .• The commutative product of 
sgn (t) g {\t\) and S (r) exists and 
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[sgn{T)g{\T\)]5{r) = 
for arbitrary g {\t\) bounded on every finite interval. 

Proof. We would like to show that 
{[sgn{T) g {\t\)]S (t) ,(p) — 0. Using the mean value 
theorem ip (t) — ip (0) + rdip (^r) /dr with < ^ < 1, we 
have 



Appendix A: Propositions on the product and the 
composition of distributions 

To investigate how the Friedmann universe crosses a 
soft singularity, we must solve the field equation in dis- 
tributional sense. For this purpose we give the defini- 
tions of the product and of the composition of distri- 
butions and prove two propositions. Fisher derived the 



following result 



sgn I 



6{r) 



for A > -1 



Our first proposition generalizes this equation for 
A > 0. The second proposition generalizes Antosik's re- 
sult: [1 + S{t)]~' = 1 [5li- Finally, we show a corollary. 

Let p (r) be any infinitely differentiable function hav- 
ing the following properties: i) p{t) =0 for |t| > 1; 
ii) p{t) > 0; Hi) pir) = p{—t); iv) j_iP(j)dT — 1. 
Then (5n (r) — np{nT) (with n = 1, 2, ...) is a reg- 
ular sequence of infinitely differentiable functions con- 
verging to Dirac delta function: lim„_^oo {^m v) = (^i v) 
for any (p £ V [s^]- Here V denotes the space of test 



< 



l(Nn(T).g(|T|)]„5„ (r),^)I 
^(0)/ [sgniT)g{\T\)]^Snir)dT 

J~l/n 

dip (r) 



dr 



sup 

|r|<l/ra 
fl/n 

X / \T[sgn{T) g {\T\)\^5^{T)\dT . 

-1/n 



The first integral on the right side of the above equation 
vanishes because the integrand is an odd function. For 



^ This definition can be generalized for the cases when the usual 
limit does not exist by taking the so-called neutrix limit [52j . 
[53| . However, we do not need for this more general definition 
here. 
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the second integrand, we have 



\T[sgn (r)5(|r|)]„(5„ {T)\dT 

-l/n 

fl/n /"l/^ 

V^nir)] j \g{\T - x\)\5n{x)dxdT 

'-l/n J-l/n 

< n sup |p(r)| / |r(5„(r)| / \g {\t - x\)\ dxdr 

|r|<l/ri J-l/j^ 



l/n 



<.i/n 

< 2 sup \p(t)\ sup |5(|r|)| / |T,5„(r)|dT 

T|<l/n |T|<l/n J-l/n 



= - sup |p(r)| sup |5(|r|)| / \yp{y)\dy 

" |r|<l/ri |-r|<l/n J-1 

2 

< - sup |p(r)| sup |5(|r|)| , 

" |T|<l/ri |-r|<l/n 

that vanishes for n — > cxd. 



□ 



Definition 2. ; T/ie composition F (/) o/ distributions 
F and f exists and is equal to h (z V on the interval 
{a,b) if 



hm 

n— ^oo 



hm / -F„(/m(r)) (p (t) dr 



{h,ip) 



for all if £ V with support contained in the interval (a, b) 

2 



Proposition 2. ; The composition of distribution 
P{t)[R{t)+Q{t)S{t)]~'^ (where w > 0, P (r) is 
bounded, R (r) 0, and in some range close t = the 
signs of R (r) and Q (r) are the same if Q (j) ^Q) exists 
if P (r) /i?"^ (r) exists^ and 



P{r) 



P{r) 



[R{T)^Q{T)5{T)r R-{r)- 

Proof. By the definition of composition of distributions, 
we should calculate 



Pnir) 



[R^ (r) + Q,n (r) (5™ (r)];; 
Performing a change of the variables as t — r, y 



^ This definition can be generalized for the cases when the usual 
limit does not exist by taking double neutrix limit |54| , |55ll , |56|| . 

^ We note that this proposition can be held even if P (t) = 1 and 
-R(t) = 5 (t) with = 1, 2, .... Indeed, 6~'^ (t) exists in neutrix 
limit and <5~" (t) = [55j . Thus in the definition 2, the usual 
limit must be changed for neutrix limit for this case. 



m{T — x), we have 



oo />oo 



1 

m 
1 

./-oo Jill 
/•oo 



-oo — oo 
oo 



1 

m 



LP (r) Pn (r) ^„ (t - y/m) 

[Rm (y/m) + mQm {y/m) p [y)] 

If (t) P„ (t) (5„ (r - y/m) 
-, — ; — ^ dydT 

R^ (y/m) 

'PjT) Pn{T)Sn {t - y / m) 

n2 (y/ + "^0™ (y/™) p (y)]' 



where = {y ■ \y\ < 1 and p (y) ^ 0} and fli 
The double limit of the first term is 



-dydr 



rdydr , 



no 



lim lim 

n— >oo m— ^— ^oo 777, 



dTcp (r) P„ (r) 



^n (t - y/m) 



dy 



R'i'niy/m) 

dT</j (r) P„ (r) 

O 

5„ {x) 



= lim lim 

n^oo m— ^— ^ 



m I r — 2; I G i i 1 



P{r) 
R'^ (r) 



We investigate the absolute value of the second integral. 
According to our assumptions for R and Q, and since we 
are interested in tti — > oo, we can choose m large enough 
to let the signs of R and Q be the same, then for a; > 0: 



< 



7I+" 



ip (t) P„ (t) Sn (r - y/m) 
02 [Pm (y/m) + mQ,n {y/m) p (y)] 

/ X p / X , Sn{T-y/m) 
LP (t) Pnir) 



n2 Qm {y/m) (y) 



-dydr 



dydr 



Performing a change of the variables as z = n (t — y/m), 
y = y, we have 



< 



< 



1 



1 



1 Jn 



sup 



Vri mJ \n m 



y\ P{z) 



ml (y) 

y 



dydz 



m- ■ - n2,|z|<i 

r-l ^1 

X / p{z)dz dy 



^(- + ^)pJ- + ^)p-iy) 

\n m/ \?7- m/ 



-1 



sup 



m- ' f22,|z|<i 
that vanishes for m 
1. 
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— I P 

n m/ 



(y) 



cxD if P is bounded. 

The 



□ 

distribution 



Corollary 

d|p(r)[P(T) + Q(r)J(r)]""|/dT (with the same 

properties for P , R, Q and oj as in proposition 2) exists 
if P (r) /P" (r) and its derivative exist, and 



P{t) 



d Pjr) 

dr [P(r) + Q(t)(5(t)]" dr R" {t) 
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Proof. Applying the derivative of a distribution at tests 
functions, and using the fact that dtp/dr £ V for any 
(y9 G 2?, and by the proposition 2, we have 



P(r) 



dr [R (r) + Q{t)6 (r)] 

Pjr) 

[R{T) + Q{T)S{T)r' dr 



d 



Pjr) 

R'^ (r) ' dr 



d P{r) 
dr R'^ (r) 



□ 



Appendix B: Two simple examples of the product 
and of the decomposition of distributions 

The definition of the product and the composition of 
distributions, used in this paper and presented in the 
Appendix A are not often encountered in physics. Thus, 
to give the reader some flavor of the corresponding con- 
siderations, using simpler means we decided to give two 
semi-heuristic examples of such products and composi- 
tions. We consider first a remarkable formula 



(Bl) 



which was first proven in Here P means the princi- 
pal value of the corresponding function. We shall prove 
here that the regularizing succession of functions with 
compact support p, employed in the Appendix A and 
the references therein, can be chosen alternatively as the 
family of the Cauchy-Lorcntz functions 



feix) 



(B2) 



It is well known that when the small parameter e — >■ 0, the 
functions of this family tend in the distributional sense 
to the Dirac 6 function. Obviously, the convolution of 
the function (jB2l with the Dirac 6 function gives again 
the same function (IB2I): 



(B3) 



The calculation of the convolution of the principal value 
P {■^) with the function fe{x) is slightly more compli- 
cated: 



P{ - ]*Mx) = lim 

X I e^O 



dy 



1 



dy- 



X — y Tr{y^ + e^) 

X 



c+e " X ~ y Tr{y^ + e"^) ^ 
The product of the expressions (|B3p and (|B4p is 



(B4) 



(B5) 



Let us now consider a family of functions 



dfejx) 

dx 



2xe 



'k{x'^ 



(B6) 



One can easily prove that if the family of functions (jB2[) 
converges in the distributional sense to the Dirac (5 func- 
tion, the family of their derivatives (|B6P converges to 
the derivative of the delta function. Now, comparing the 
right-hand sides of Eqs. (jB6l) and (jB5|) we see that when 
e — > the product in the left-hand side of Eq. (jB5|) con- 
verges in the distributional sense to —\S'[x) and thus the 
correctness of the equality (jBl[) is checked. 
Now let us discuss the Antosik identity j5l| 



1 



l + (5(a;) 



1. 



(B7) 



Here we have the composition of the distributions ^(5), 
where F = ^ and (/(x) = 1 -I- i5(a;). Calculating the con- 
volutions of the distributions F and g with the Cauchy- 
Lorentz functions (IB2D we obtain 



,{x) = 1 + 



5 



(B8) 



(B9) 



Correspondingly the composition of these functions is 

1 



Faig.) 



^2+1 



and it is easy to check that 



lim lim F^fg^) = 1, 

CT^O e-i-0 



confirming the identity (|B7 



(BIO) 



(Bll) 



Appendix C: The Lanczos equation 

For a generic junction surface the Lanczos equation 
emerges from the Gauss-Codazzi relations [HI], [i^. The 
projected Lie derivative of the extrinsic curvature Kab in 
the normal direction n to the surface is 

KhiCnK,, = -3e (^KhtT^k - ^g'^T,,^ + Zab (CI) 

(Eq. (21) of dl] in the units SttG/S = 1), with 

2^ab — -eP-ab + "^KacK^j ~ g^^KikKab 

+Dbaa - eabUa ■ (C2) 

Here gab is the space-time metric, hab ~ gab ^ ^^aUb 
(e = n'^Ua — {^1; 1}) is the induced metric on the junc- 
tion surface, and Tab is the energy-momentum tensor. 
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The tensor Zab depends only on geometrical quantities: 
TZab and D are the Ricci tensor and covariant derivative 
induced on the hypersurface, and Ua — "nf^Vcna, with V 
the 4-dimensional covariant derivative. 

When the energy-momentum tensor is a sum Tik — 
Hik + TikS (t) (where t is the coordinate adapted to n, 
i.e. n = tg^d/dr, and n°-Tab = 0), with 11^^ the regular 
4-dimensional part and T^^ the distributional part on the 
hypersurface, integration of Eq. (jCl|) across r through 
an infinitesimal range containing the hypersurface keeps 
only the distributional part, leading to the Lanczos equa- 
tion (13, |40i. 



^Kab = -3e ( Tab - -^^ab 



or equivalently 



3eTab - AKab - habAK 



(C3) 



(C4) 



Let us now specialize this for a junction along a maxi- 
mally symmetric r = spacelike hypersurface (a hyper- 
plane with TZab = 0) embedded in a flat Friedmann space- 
time. The normal vector n has zero acceleration = 
and the extrinsic curvature becomes Kab — aahab, with 
hab the 3-dimensional^ Euclidean metric. The curvature 
term is Zab = —H'^a^hab and the energy momentum ten- 
sors are Hab — 'p^aUb +pa?hab and Tab = pa^hab- Since 
the projected Lie-derivative in Eq. (jCip becomes a time 
derivative, the equation reads 



(C5) 



which is a combination of the Raychaudhuri and Fried- 
mann equations. For finite H, p and p as before the in- 
tegration of Eq. (jC5l) across an infinitesimal time range 
T leads to the Lanczos equation 



Here T is the trace of Tab- As Zab is finite, its contribu- 
tion to the integral across the infinitesimal range also van- 
ishes. Without a distributional energy-momentum part, 
the extrinsic curvature should be continuous. 



(C6) 
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